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Abstract 

Three years after the completion of the next-to-leading order calculation, the status of the 
theoretical estimates of e'/e is reviewed. In spite of the theoretical progress, the prediction of e'/e 
is still affected by a 100% theoretical error. In this paper the different sources of uncertainty are 
critically analysed and an updated estimate of e'/e is presented. Some theoretical implications 
of a value of e'/e definitely larger than 10 -3 are also discussed. 
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1 Introduction 


After many years, the direct CP violation in K 0 decays, parametrized by e', is still an open 
issue. The last generation experiments have found [l], [2|] 

£ f (7.4 ± 5.9) x 1CT 4 E731 m 

e “ l (20 ± 7) x 10" 4 NA31 ^ ’ 

From these results no definite conclusions can be drawn on the CP property of the K° decay 
vertices, namely on whether e' is vanishing. 

Theoretically, the Standard Model makes precise assumptions on the mechanism that gener¬ 
ates the CP violation. Indeed the only source of CP violation is the free phase which appears in 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix [|3jj with three quark generations. This choice 
implies a non-vanishing e', unless some dynamical cancellation occurs. The CP-violating phase 
appears in the K° decays through the so-called penguin diagrams, ffowever other choices are 
possible, for example the superweak model 0, which predicts strictly e' = 0. 

To clarify this issue, a new generation of experiments is going to be built, achieving a 
sensitivity on e'/e at the level of 1-2 x 1CU 4 ||. On the theoretical side, the problem is giving 
a reliable estimate of e'/e, including the theoretical error. The task is not easy: physics from 
many scales effectively contribute to e'/e, from the top mass down to the strange mass, including 
important non-perturbative effects. Nevertheless all recent analyses agree on predicting e'/e = 
few x ICO 4 , with roughly a 100% relative error Sla¬ 
in the following we present an updated prediction of e'/e, giving an account of the procedure 
and the different sources of theoretical uncertainty. We also discuss the dependence of e'/e on 
some critical non-perturbative parameters. 

2 e'/e in a few steps 

The essential theoretical tool for the calculation of e'/e is the AS = 1 effective Hamiltonian, 
which allows the separation of the short- and long-distance physics. Using the effective Hamil¬ 
tonian, one obtains an expression of e'/e that involves CKM parameters, Wilson coefficients and 
local operator matrix elements. Therefore the evaluation of e'/e requires essentially three steps, 
namely (1) the phenomenological determination of the CKM parameters, (2) the calculation 
of the Wilson coefficients at a next-to-leading order (NLO) and (3) the determination of the 
matrix elements of the local operators appearing in the AS = 1 effective Hamiltonian. 
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2.1 Basic formulae 


The NLO AS = 1 effective Hamiltonian at a scale mj, > /i > m c can be written as 
A,.Gj 


n = 


7T 


{(1 - r)^)^) - Q5(/x)) + C 2 (/i) (Q a (/i) - QSM) 


+^Ec ! (/i)g i (/i)} , 


( 2 ) 


i=l 


where is the Fermi constant, A g = V q dV* s and r = —X t /X u , V qiqj being the CKM matrix 
elements. The CP-conserving and CP-violating contributions are easily separated, the latter 
being proportional to r. 

The operator basis includes eleven dimension-six local four-fermion operator^]. They are 
given by 


Q 1 = {s a d a )(v-A){upUp){y-A) 

Q 2 = (s a dp)[y-A){upU a )iy-A) 

Q 1 { s ada) (V—A){ c (3Cp) (y-A) 

Q C 2 = ( S a dp)y-A)(cpC a )(y-A) 

Q3,5 = (Sqdg) (V—A) 

9 

< 34,6 = ( S a dp)(V-A) y^XQi3Qa)(VTA) ( 3 ) 

9 

3 

^ 7,9 = 7 ;(Sgda)(V-A) ^aiQgQa) (V±A) 

Z 9 

3 

^8 = -(s Q d/3)(y-A) e q {qpqa){V+A) , 

Z 9 

where ( q a Qp)(v±A) — <?«7m(1 4 = 75)^, a and /3 are colour indices, and the sum index q runs over 
{d, u,s,c}. Operators Q3-Q5 are generated by the insertion of the tree level operator Q 2 into 
the strong penguin diagram, while Q7-Q9 come from the electromagnetic penguin diagrams. 
As we will see, the two classes of operators are both relevant for e'/e. Further details on the 
NLO AS = 1 effective Hamiltonian can be found in ref. || 

From the definition of e' and using eq. (^]), one readily obtains 


/ .e 
e = 1 


i(S2—So) 


CO 


y/2 ReA 0 


a; 1 (ImA 2 ) / — (1 — O/s) Iim4 0 


1 One more operator must be included if /i > mb is considered. 


( 4 ) 
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where, as usual, A/e lSl = {'K'k{I)\H\K g ) , Qjb is the isospin breaking contribution due to the 
Tr-rj-rj' mixing ||, uj — ReA 2 /ReAo, and 


ImAn = 


-G F i™(y; s v u ){-(c & B & + l -c b Bl)z + (c a b a + 1c 3 b 3 )x 
+C 7 By\^ + f + f) + C 8 Bl'\2Y + f + f) 


-C 9 B 9 


1/2 


f + (^ + c 2 fl S )x}, 


(5) 


(ImA 2 )' = 


-a F im(y t yu) {c 7 b‘ /2 (h - A) + c 8 b 8 3/2 (k - A 

+c sB rf} • 


The relevant operator matrix elements are given in terms of the R-parameters as follows: 

(7nr(0)\Qi\K) = Bl 2 (Tnr(0)\Qi\K) VIA 

(7rn(2)\Q i \K) = B* /2 (nir(2)\Q t \K) VIA , (6) 


where the subscript VIA means that the matrix elements are calculated in the vacuum insertion 
approximation. VIA matrix elements can be calculated and expressed in terms of the three 
quantities 


X 

Y 

Z 



fn 


M 2 K 


.m s (n) +m d (n ) 



/150 MeV\ 

V m s{lA ) 


( 7 ) 


Notice that, contrary to X and Z , Y does not vanish in the chiral limit. This reflects the 
different chiral properties of the operators Q 7 and Q 8 . 


2.2 CKM matrix elements 


In order to estimate e'/e, we need lm\V^ s V t( ij. This requires a complete knowledge of the CKM 
matrix V. For example, in the Wolfcnstein parametrization [0, up to 0(A 3 ), 


V 


( 1 — 4r A AA 3 (p — ir]) 

-A 1 - f AA 2 

\ AA 3 (1 — p + irj) —AX 2 1 


( 8 ) 
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m t = 1 67 ± 3 GeV A W QCD =330± 1 OO MeV 



77 vs p rj vs p, f B cut 


Figure 1: Constraints on the CP-violating phase and the unitarity triangle. The cos 5 distribu¬ 
tions are plotted above. The dashed distribution includes the constraint coming from AMB d and 
the lattice determination of Bsy/Jb- Below, the same constraints are shown as contour plots 
in the p^rj plain. The solid, dashed and dotted lines correspond to the 5%, 68% and 95% of 
the generated configurations respectively. In this plane the allowed region determines the third 
vertex of the unitarity triangle J2 q ={ u ,c,t} VqiXqd = 0, the others being (0, 0) and (1, 0). 


and 

lm(v*V td ) = -A 2 X 5 g = -A 2 X 5 asin8 , (9) 

using the definition of the C P-violating phase 5 given by cre l5 = p + irj. While A is well known, 
A and a can be extracted from the measurements of the B lifetime and the semileptonic decay 
rates, respectively. The remaining task is the determination of the C P-violating phase. The 
CP-violating parameter e in the K°-K° mixing is the obvious tool. It is given by 

|e|g =0 = C e BxA 2 X 6 asin8 ^F(x c ,x t ) + F(x t )[A 2 A 4 (l — acostf)] — P(x c )| , (10) 
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g f 

1.16639 x 10 -5 GeV -2 

m c 

1.5 GeV 

m b 

4.5 GeV 

M\y 

80.32 GeV 

M b 

5.279 GeV 

M k 

498 MeV 

A M k 

3.495 x 10- 12 MeV 

M n 

140 MeV 

u 

132 MeV 

Ik 

160 MeV 

ReH 0 

2.7 x 10 -7 GeV 

a; 

0.045 

A = sin 9 C 

0.221 

^exp 

2.284 x 10“ 3 

d 

2 GeV 


Table 1: Input parameters assumed to be constants in the analysis 


where F(xi ) and F(xi,Xj ) are the Inami-Lim functions [ 111, including the QCD corrections 
I2[ , and 

G\flM K Ml 


c f = 


( 11 ) 


6a/2tt 2 A M k 

The comparison of the previous expressions with the measured values of e allows the extraction 
of cos 5. However, since eq. © is a quadratic function of cos 5, one obtains two different 
solutions, corresponding roughly to cos 5 positive and negative. The solid-line distribution of 
cos S in fig. |l] is the result of this analysis. The distribution is obtained by generating randomly 
the various parameters according to the values and errors listed in tables [3] to |3]. 

To further constrain the GP-violating phase, one can exploit the Bd — Bd mass difference, 
given by 

A M Bd = {^B b A 2 A 6 (l + a 2 - 2cr cos S^jF(x t ) , (12) 

where B b is the P-parameter associated to the A B = 2 operator (bd)(y_A){bd)(y_A)- Recent 
lattice results give quite large values of PbVTb EH- We use Bb\[Jb = 210 ± 35 MeV. The 
requirement of compatibility between B b ^/Jb extracted from eq. ([12]) and the lattice value 
results in an effective selection of the positive values of cos <5; see the dashed cos 5 distribution 
in fig. |l]. Therefore we can quote the value 


cos 5 = 0.38 ±0.23 , 


(13) 
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1V QCD 

330 ± 100 MeV 

fllB 

0.25 ±0.10 


210 ± 35 MeV 

ryRG—inv 

b k 

BU 

-£>3,4 
-S.5,6 

-°7-8-9 

b \ 3/21 

j^(3/2) 

o(3/2) 

n 9 

0.75 ±0.15 

0 - 0.15 w 

1 - (>M 

1.0 ± 0.2 

1 « 

0.6 ± 0.1 

0.8 ±0.15 

0.62 ± 0.10 


Table 2: Input parameters assumed to be uniform in the analysis 
where the error is the variance of the dashed distribution in fig. |I]. 

2.3 Wilson coefficients 

The Wilson coefficients Cfp) appearing in eq. (|| can be calculated using the renormalization 
group improved perturbation theory, provided that p is a scale large enough for the perturbation 
theory to be reliable. Indeed their p dependence is controlled by the renormalization group 
equation 

ttfhc i ( P ) = l -Y.%i c M, (u) 

where 7 is the anomalous dimension matrix of the operators in eq. (£IJ). The NLO calculation 
of these coefficients is discussed in ref. [[|. To our end, it suffices to recall that, for any suitable 
value of / 1 , the Wilson coefficients are a known set of real numbers, which however still depend 
on the scheme chosen to renormalize the local operators. The typical relative error on the 
coefficients relevant for e'/e is 10-20%. For more details and numerical values, see for example 
ref. [§. 

2.4 Local operator matrix elements 

The calculation of the matrix elements of the local operators appearing in eq. (2j) requires 
the use of a non-perturbative technique. Indeed these matrix elements contain the low-energy 
QCD dynamics, from the scale p downward. Besides the p dependence, at the NLO they also 
depend on the operator regularization scheme. Both these dependences must be matched with 
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167 ± 8 GeV 

\V c b\ = AA 2 

0.040 ± 0.003 

\V ub /V cb \ = Xa 

0.080 ±0.015 

tb 

1.56 ± 0.06 ps 

A M e B x f 

0.464 ± 0.018 ps" 1 

mf( 2 GeV) 

128 ± 18 MeV 


Table 3: Input parameters assumed to be Gaussian in the analysis 


the corresponding dependence in the Wilson coefficients, in order to have a scale- and scheme- 
independent physical prediction. Therefore only a non-perturbative approach which allows a 
full control over the renormalization scale and scheme dependences can be consistently used. 
Furthermore this technique should allow choosing the renormalization scale fi large enough 
for the perturbative calculation of the Wilson coefficients to be reliable. As far as e'/e is 
concerned, the only known non-perturbative approach that fulfills these requirements is lattice 
QCD. Several 5-parameters appearing in eq. (||) have been computed on the lattice [|14|]. In 
particular 5 6 and 5g 7 , which turn out to be the most important numerically, are known. 
For the others, we use the vacuum insertion approximation, namely B — 1 at fi = 2 GeV. 
There are two exceptions: first, 5 3 and 5 4 are allowed to be as large as 6, considering that the 
penguin operator matrix elements may be at least partially responsible for the A I = 1/2 rule. 
The other is 5 42 , which strictly speaking cannot even be defined, since (7r7r|Q 4 2 |Ji°)y/^ = 0. 
However, a small contribution is expected beyond the VIA and we parametrize it by assuming 
the VIA matrix elements to be equal to those of Qi ,2 and introducing a small 5-parameter. 
Table [2] contains the values of the 5-parameters used in the numerical analysis. 

The relevant 5-parameters are known with an error of about 20%. We will see that unfor¬ 
tunately they produce a larger error in e'/e because of a partial cancellation between different 
terms. However, there is another source of uncertainty stemming from the normalization of the 
matrix elements, namely the value of the running strange mass, see eqs. The former 

good agreement between lattice and QCD sum rules on this mass 0 |TT)|| has been questioned 
by a recent lattice result, in which the extrapolation to the continuum limit is attempted [ P~8[| . 
We will quantitatively discuss this issue in the next section; however, it is worth noting that 
this uncertainty comes only from the choice of normalizing the lattice results to the vacuum 
insertion approximation. In the future the possibility of normalizing lattice results to a better 
known quantity should be considered. 
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m t = 167±8 GeV 
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Figure 2: Plot of the distribution of e'/e, p(eVe), above, and of Pie'/e > x) 
JT dyp(y)/ S- OQ dyp(y), below 


2.5 The “best” estimate 

Having the necessary ingredients, we can put them togheter and produce an estimate of e'/e. 
Varying randomly the parameters in tables |2] and |3], we obtain the distribution of fig. from 
which we estimate 


= (4.6 ±3.0 ±0.4) x 10 


-4 


(15) 


Again the first error is the variance of the distribution, while the second one refers to the 
residual scheme dependence coming from higher order in perturbation theory and it is obtained 
by using two different renormalization schemes. 

Figure H also shows the probability distribution P(e'/e > x) as a function of x. For instance, 
given our choice of the input parameters and distributions, one obtains P(e'/e > 2x 10 - 4 ) ~ 0.8. 

Why is the theoretical error 100% or even larger? Wilson coefficients and P-parameters 
have errors ~ 20% or less. One can argue that there are many contributions, but fig. fl shows 
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c 1/2 3/2 

Figure 3: Relative contributions of the different terms appearing in eq. (W with respect to e'/e. 

(a /o') 

The largest terms are those containing B$ and B$ ' ’. A partial cancellation occurs between 
them. 

that actually the bulk of the result can be obtained retaining only two operators. On the other 
hand the two main terms, say Bq and B s , partially cancel each other, lowering the central 
value of the prediction and increasing the relative error. The effectiveness of this cancellation 
depends not only on the top mass, which nowadays is a well-known quantity, but also on 
the isospin-breaking parameter fljB, of which we have only quite old theoretical estimates |J. 
Contributions to the error also come from some overall factors, sin 8 and mf 2 which appears in 
front of the largest terms. All these effects sum up to give the large error of the final result. 

3 The future of e'/e 

Given the theoretical difficulties in reducing the error on e'/e, we may wonder what we will learn 
from the next generation of experiments, with an expected sensitivity at the level of 2 x 10~ 4 . 

3.1 What the present estimate is good for 

According to present theoretical estimates, the Standard Model can accomodate e'/e values 
between 0 and 10 -3 . If the present ambiguity on the experimental value of e'/e will be solved, 
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confirming the E731 figure, then the new measurement will hardly contribute to improving our 
knowledge of the Standard Model and the QCD dynamics. Indeed, even if a non-vanishing value 
of e'/e would rule out the superweak model allowed at present, a measurement of e'/e < 1CT 3 
can be accommodated by the Standard Model with no pain, weakly constraining the different 
parameters that enter the theoretical expression of e'/e. 

On the contrary a measurement of e'/e > 10~ 3 would be a fruitful surprise. Its explanation 
would require either new physics at work or an improved understanding of the QCD dynamics. 


3.2 What does e'/e > 10 3 imply? 


In the following let us assume that indeed the future measurements of e'/e give a value definitely 
larger than 10~ 3 . New sources of CP violation beyond the Standard Model can change the 
prediction in eq. (1TB) and accommodate such a large value. However the task is not trivial. 
For instance supersymmetric effects on e'/e are small in comparison with the large error, in the 
Minimal Supersymmetric Standard Model |jT7l| . 

Anyway, let us ignore here the possibility of new physics effects and stick to the Standard 
Model. One or more input parameters must deviate from the values listed in tables M to 
account for the increase of e'/e in this case. We believe that these tables collect the “best” 
set of input parameters allowed by our present knowledge of the CKM matrix and the QCD 
dynamics. However, one can still consider the effect on e'/e caused by assuming that some 
critical parameter is outside its allowed range. This is a particularly reasonable speculation in 
the case of the H-parameters and the hadronic quantities in general. Indeed the systematic 
errors of the lattice QCD results are not completely under control at present and may be 
underestimated in some cases. Hopefully new developments of the lattice QCD techniques, 
such as non-perturbative renormalization and unquenching, will lead to a better control of the 
systematic errors in the future. 

Let us consider two scenarios leading to a large e'/e: first a small strange mass, then 
anomalous B$ and B 8 . 

The possibility of a small strange mass, at the level of 50-90 MeV, has been suggested by 
a recent compilation of lattice results |T|. A previous lattice analysis, done at finite lattice 
spacing, gave a running mass mf rs (2 GeV) = 128 ± 18 MeV [15], in agreement with the 
QCD sum rules determination |16|. The new analysis of ref. |18j extrapolates the value of the 
strange mass to zero lattice spacing, finding m / rs (2 GeV) = 90 ± 15 MeV in the quenched 
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approximation. The inclusion of dynamical fermions further decreases m s and they quote 
m,g IS (2 GeV) = 70 ± 11 MeV for two generations of dynamical fermions. This result indicates 
at least that the error associated to the finite lattice spacing was probably underestimated in 
the previous analysis. However, we believe that it is premature to accept the new figures for 
the running strange mass, until the discrepancy with the QCD sum rules result is understood. 

Anyway it is an easy exercise to see how e'/e changes in this small-m. s scenario, since the 
largest terms are proportional to m~ 2 . We obtain 


- = (14 ± 8) x 10 
e 


-4 


(16) 


for (2 GeV) = 70 ± 11 MeV. Thus a small strange mass can push e'/e up to some units in 
ICE 3 . Notice however that the relative error is roughly unchanged and a vanishing e'/e is still 
well inside the allowed range. We stress again that, as far as e'/e is concerned, the problem of 
m s appears as the consequence of an inappropriate choice of the H-parameter normalization. 

In another conceivable scenario, one can take Bq to be large and at the same time a small 
B^ 2 . This choice tends to spoil the cancellation between strong- and electro-penguin operators, 
thus increasing e'/e. Indeed a recent result gives a new value B^ 2 = 0.81 ± O.OSiJ/o 3 smaller 


than the previous ones |]18|, [Hj. To our knowledge, no new results are available for B e . We can 
explore an extreme scenario where we choose 


B„ = 1.50 ± 0.15 and B% /2 = 0.50 ± 0.05 . 


(17) 


In this case, we obtain 

- = (12 ±4) x 10“ 4 (18) 

Again e'/e larger than ICE 3 is predicted. The relative error is smaller than before because of the 
spoiling of the cancellation between B 6 and B s . Now e' = 0 is only marginally compatible 
with the prediction. Notice, however, that precisely the very peculiar situation of a large B$ 
and a small B g 3 ^ 2) must be realized in order to enhance e'/e. 

The previous examples show that a measurement of e'/e larger than ICE 3 would indeed have 
non-trivial implications, forcing us either to reconsider some non-perturbative results or to call 
for new physics. 
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